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ABSTRACT 

We study the behaviors of galactic disks in triaxial halos both numerically and 
analytically to see if warps can be excited and sustained in triaxial potentials. We 
consider the following two scenarios: 1) galactic disks that are initially tilted relative 
to the equatorial plane of the halo (for a pedagogical purpose), and 2) tilted infall of 
dark matter relative to the equatorial plane of the disk and the halo. With numerical 
simulations of 100,000 disk particles in a fixed halo potential, we find that in triaxial 
halos, warps can be excited and sustained just as in spherical or axisymmetric halos 
but they show some oscillatory behaviors and even can be transformed to a polar-ring 
system if the halo has a prolate-like triaxiality. The non-axisymmetric component of 
the halo causes the disk to nutate, and the differential nutation between the inner and 
outer parts of the disk generally makes the magnitude of the warp slightly diminish 
and fluctuate. We also find that warps are relatively weaker in oblate and oblate-like 
triaxial halos, and since these halos are the halo configurations of disk galaxies inferred 
by cosmological simulations, our results are consistent with the fact that most of the 
observed warps are quite weak. We derive approximate formulae for the torques exerted 
on the disk by the triaxial halo and the dark matter torus, and with these formulae we 
successfully describe the behaviors of the disks in our simulations. The techniques 
used in deriving these formulae could be applied for realistic halos with more complex 
structures. 

Subject headings: galaxies: evolution — galaxies: halos — galaxies: kinematics and 
dynamics — methods: n-body simulations 



1. INTRODUCTION 

The outer parts of the disks in most spiral galaxies, even in isolated ones, are warped (Ann 
& Park 2006 and references therein), and this implies that the warps must be either long-lived 
or be repeatedly excited. Early proposals for the excitation mechanism of disk warps include the 



x Dept. of Astronomy & Space Science, Kyung Hee University, Yongin-shi, Kyungki-do 449-701, Korea; my- 
jeon@ap4.khu.ac.kr, sungsoo. kim@khu.ac.kr 

2 Division of Science Education, Pusan National University, Pusan 609-735, Korea; hbann@pusan.ac.kr 

3 Corresponding author 



- 2 - 



intergalactic gas flow past the disk (Kahn & Woltjer 1959), tidal distortion due to the Magellanic 
Clouds (Elwert & Hablick 1965), and the free mode of oscillation of a disk (Lynden-Bell 1965). 

Extending the work by Lynden-Bell, Hunter and Toomre (1969) studied the warp in terms of 
discrete modes of oscillation of an isolated, thin, self-gravitating disk and found that a long-lived 
warp can correspond to a simple discrete mode only when the disk has an unrealistically sharp 
edge. Later, Sparke and Casertano (1988) showed that self-gravitating disks even with realistic 
density profiles may have long-lived discrete modes if the disk is embedded in the fixed potential 
of an axisymmetric halo (such modes are called the modified-tilt mode). 

Dekel & Shlosman (1983) and Toomre (1983) suggested a misalignment between the disk and 
the non-spherical halo as the cause of the warp, and Dubinski & Kuijken (1995) and Ideta et al. 
(2000), among others, performed simulations for such configuration. With simulations of tilted 
disks in fixed, axisymmetric halos, Ideta et al. found that prolate halos could sustain galactic 
warps while the warping in oblate halos continues to wind up and disappears fairly quickly. 

However, dynamical friction between the disk and halo may cause the warp to disperse within 
timescales much shorter than a Hubble time (Nelson and Tremaine 1995). Numerical studies em- 
ploying live halos (Dubinski & Kuijken 1995; Binney, Jiang, & Dutta 1998), which can properly 
consider the reaction of the halo to the disk, confirmed that dynamical friction may play an impor- 
tant role in the disappearance of galactic warps on a time scale shorter than a Hubble time. Thus 
a successful scenario for the warp should not give rise to significant dynamical friction. Plausible 
scenarios that do not suffer the damping problem include the reorientation of a massive galactic 
halo by the cosmic infall (Ostriker & Binney 1989; Jiang & Binney 1999; Shen & Sellwood 2006) 
and the external torque produced in the disk by the accretion of the intergalactic medium (Revaz 
& Pfenninger 2001; Lopez-Corredoira et al. 2002). 

Shen & Sellwood (2006) performed simulations of an idealized form of cosmic infall on to a 
disk galaxy and obtained warps that persist for a few Gyr. They found that the damping from 
the live halo is very weak compared to that in initially tilted disk scenarios, and attributed the 
difference to the facts that 1) the source of the torque in their model resides in the outer halo only, 
2) the precession rate of their disk is very low, and 3) the effect of the torque in their model kicks 
in gradually. 

It appears that the model by Shen &; Sellwood (2006) is a highly feasible way of exciting and 
sustaining a warp, because their warps closely resemble those observed and do not significantly 
suffer the damping problem. But their halo is modeled to be spherical, wherease many numerical 
studies on the shape of the dark matter halos anticipate that the triaxiality of the dark matter 
halos is not negligible (Jing & Suto 2002; Bailin & Steinmetz 2005; Allgood et al. 2007, among 
others). Thus it is worth checking if the results of the cosmic infall scenario are significantly altered 
when the halo is modeled to be axisymmetric or triaxial. Furthermore, as the triaxiality of the 
halo exerts a torque that periodically increases and decreases the inclination angle of the tilted disk 
(see §2j), the initially tilted disk scenario deserves a revisit to see if the extra torque by the halo 
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triaxiality can diminish the effect of the damping. 

For these reasons, in the present paper we study the evolution of the self-gravitating disks in 
axisymmetric and triaxial halos for two scenarios: 1) initially tilted disks in triaxial halos and 2) 
tilted cosmic infall in axisymmetric and triaxial halos. As discussed above, the first scenario does 
not appear to be a plausible mechanism for real warps, but we first study this scenario in detail to 
understand the influence of the triaxial potential alone to the disk. After successfully describing it 
with simulations and some analytical formulae, we then analyze the behaviour of the warp in the 
second scenario. 

As a first step, the halos in the present paper are fixed, and this will enable us to concentrate 
on studying how the axisymmetry and/or triaxiaility of the halo alters the evolution of the warp 
compared to the spherical or axisymmetric cases. As a subsequent study to the present one, we 
will examine the disks in "live" axisymmetric and triaxial halos, which can properly treat the effect 
of dynamical friction between the disk and the halo. 

Note that the aim of the present study is to understand the fundamental dynamics between 
the disk and the axisymmetric or triaxial halo for the above two scenarios, rather than to find the 
conditions that best reproduce the various observed characteristics of the warp. Still, we will briefly 
compare our simulations and observed warps in §61 

Warps are generally more extended in radio observations of gas than in optical observations of 
stars. However, in the present study, we only consider the stellar component of the disk to avoid 
the model complexity and computing cost involved with the consideration of the hydrodynamic 
(and even magnetohydrodynamic) effects. Thus the warps in our simulations are to be compared 
to the inner parts of the observed, extended warps. 

We describe the models and the numerical method in §21 and discuss some theoretical back- 
grounds necessary for the analyses of our simulations in §3l The results and analyses of our sim- 
ulations for initially tilted disks are presented in §H and those for the tilted cosmic infall scenario 
are presented in §5l We summarize and discuss our results in §61 

2. MODELS 

For the simulations presented here, we use a parallel version of the tree N-body code named 
GADGET (Springel, Yoshida, and White 2001). The code can calculate hydrodynamic forces using 
the smoothed particle hydrodynamics (SPH) technique, but we only consider gravitational forces 
for the present work, i.e. we do not consider the effects of gas component in the disk. Gadget was 
the choice of numerical method also for a recent study on tidal structures in a disk galaxy created 
by gravitational interactions with a perturbing companion (Oh et al. 2008). 

Our model galaxies consist of a disk and a dark matter halo. The disk is represented by 100,000 
particles, and the halo is modeled by a fixed potential. Thus the evolution of the disk is governed 
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by the self-gravity within the disk and the external gravity from the halo (plus the dark matter 
torus for the cosmic infall scenario), but the effect of the disk on the halo and its reaction back to 
the disk, such as dynamical friction, is not considered. In this way, we would be able to observe 
the role of various triaxial halos on the excitation and maintenance of the warp separately from the 
effect of mutual back reaction. 

Following Ideta et al. (2000), we adopt the compound galaxy model by Hernquist (1993) and 
modify it for triaxial halos. The distribution of the disk particles is given by 

M d 2 ( % 



p d (R, z) = exp(-R/R d )sech 2 — , (1) 

4nRp d \z d J 

where R is the galactocentric radius projected on to the galactic plane, z the vertical height from 
the plane, M d the disk mass, R d the radial scale length, and z d the vertical scale thickness. The 
density profile for the triaxial halo is given by 

where is the halo mass and [i is defined by 

2 x 2 y 2 z 2 
M=-F + ^ + ? - (3) 

The triaxiality parameters a, 6, and c are the scale lengths along the x, y, and z axes. Our 
triaxial halos are named either "oblate-like" or "prolate-like": Oblate-like halos have major and 
intermediate axes in the equatorial plane, while prolate-like halos have minor and intermediate axes 
in the equatorial plane. Table [T] lists a, 6, and c values of our models. 

Using the method of ellipsoidal shells (Kellogg 1953; Chandrasekhar 1969; Binney & Tremaine 
2008), the halo potential can be written as 

GM h r°° du 

where 222 

a z + u b z + u c z + u 
and we calculate the force exerted by the halo by numerically differentiating this potential0 

Our simulations are in units of G = 1, R d = 1, and M d = 1, where G is the gravitational 
constant. If these units are scaled to physical values appropriate for the Milky Way, i.e. R d — 
3.5 kpc and M d — 5.6 x 1O 1O M , unit time and velocity become 1.31 x 10 7 yr and 262 km s _1 . 



1 After performing all of our simulations, we realized that Merritt & Fridman (1996) present an integral form of 
the force for our halo potential. They give the triaxial generalization of Dehnen (1993) models (both potentials and 
forces), and our halo model is a 7 = 1 case of the Dehnen model family. 
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We set Zd = 0.2 and = 17.8 for all of our simulations. These numbers result in the minimum 
Toomre Q parameter of about 1.5, thus the disk is stable against the bar instability. The orbital 
period at the half-mass radius, R ~ 1.7, is 13.4 in our system of units. We perform our simulations 
until T = 600, which is ~ 45 orbital periods at the half-mass radius of the disk. To achieve 
an equilibrium configuration, we sampled velocities from Gaussian distributions with means and 
dispersions derived from the Jeans equation following the prescription by Hernquist (1993). The 
total energy is conserved to better than 0.2 per cent. 

In our first set of simulations (initially tilted disk models), the disks are initially tilted by 30° 
with respect to the equatorial (x-y) plane of the halo in order to initiate the vertical oscillation of 
the disk particles. 

In the second set of simulations, the disks are initially in the equatorial plane of the halo with 
no systematic vertical motions, and following Shen & Sellwood (2006), we assume that the late 
infallng material in hierarchical galaxy formation models forms a gradually growing torus of radius 
R t = 10.7 Rd- Its mass increases linearly from at T = to M t = 2.5M^ at T = 200, and its 
cross-section has a Gaussian density profile with a standard deviation of 0.2R t . The torus plane is 
inclined at 15° relative to the equatorial plane of the halo, and its the ascending node aligns with 
the negative y-axis. The torus is represented by 5,000 particles that initially have circular motions 
with velocities equal to the local circular speed. We assume that the torus does not precess as 
its precession rate due to the interaction with the disk or halo is much longer than our simulation 
period and also assume that the torus is dynamically stable for the sake of simplicity. To avoid 
the instability of a cold, self-gravitating torus, we keep the torus particles at their initial positions 
during the whole simulation period. Note that our models for the cosmic infall scenario are different 
from Shen & Sellwood (2006) in that their halo density profile is based on the King model and 
their torus is located at a farther distance compared to the size of the disk (R t = 15Rd). The latter 
implies that the self-gravity within the disk will be relatively less important in our models. 



3. THEORETICAL BASIS 

Following Ideta et al. (2000), we analyze the behavior of the disk in terms of the relative 
differences in the inclination and in the line of ascending node (LON) between the inner and outer 
parts of the disk. In order for the disk to be observed as a warp, the inclination of the outer disk 
needs to be different enough from that of the inner disk, but the LONs of both inner and outer 
disks should not wind up and need to nearly coincide. 

The disk precesses when the torque applied to the disk has a non-zero LON component, and 
the precession makes the longitude of the LON, A, evolves. The inclination of the disk, z, evolves 
when the torque has a non-zero component that is normal to both the LON and the z-axis. Since 
the components of the forementioned torques are either parallel or normal to the LON, we adopt 
a rotating frame where the L-axis aligns with the LON of the disk and the P-axis is normal to 
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the L-axis in the halo equatorial plane such that when the positive L-axis aligns with the positive 
x-axis, the positive P-axis aligns with the positive y-axis (see Fig. [TJ. In this frame, the LON 
evolves due to the L-component of the torque, Tl (see eq. |A17| ), while the inclination evolves due 
to the P-component of the torque, Tp (see eq. |A16j ). 

Our disks have angular momenta whose z-component values are positive, i.e. our disks rotate 
counterclockwise when viewed from the positive z-axis (the V c values in eqs. |A16| and |A17| are 
defined to be positive for such rotations). Thus a positive (negative) Tl will induce a prograde 
(retrograde) precession, while a positive (negative) Tp will result in the decrease (increase) of the 
inclination. 



3.1. Torques by the Triaxial Halo 

Approximate formulae for the torques by triaxial halos are derived in Appendix [3] The 
characteristics of the halos are described with the coefficients J20 and J22 (see eqs. |A5j and |A6j 
for their definitions, and Table [2] for their values at r — 2 and 4; r is a 3-dimensional galactocentric 
radius), which are related to the spherical harmonic functions of degree I — 2 and order m = 0,2. 
Oblate halos have negative J20 values and prolate halos have positive J20. While both oblate and 
prolate halos have J22 = 0, triaxial halos have a non-zero J22. Thus equations (IA14I) and (IA15I) 
imply that axisymmetric halos exert Tl only, whereas triaxial halos exert both Tl and Tp. As a 
result, axisymmetric halos can cause changes in A only, but triaxial halos can cause changes in both 
i and A (see eqs. |A16j and |A17j ). 

Another consequence of the triaxial halo is the oscillations in the magnitudes of the torques. 
Equations (IA14I) and (IA15I) show that the contribution of the triaxial component J22 to both Tl and 
Tp has a sinusoidal form (cos2A or sin2A). In this subsection, we discuss the behavior of the disk 
in triaxial halos with a > b (all of our triaxial halos except those in models CI-OL7b and CI-PL7b 
have a > 6), which have positive J22 values. In oblate-like halos, the LON undergoes a retrograde 
precession (because of the negative J20), thus the absolute value of Tl decreases (increases) when 
A is in the first and third (second and fourth) quadrants. In prolate-like halos, the LON advances 
forward (because of the positive J20) and thus as in the oblate-like halos, the absolute value of Tl 
decreases (increases) when A is in the first and third (second and fourth) quadrants. These effects 
cause the precession rate to oscillate twice during one full precession of the LON {\d\/dt\ decreases 
[increases] when A is in the first and third [second and fourth] quadrants). 

As for Tp, the sign of its value is more important than the decrease or increase of the value. In 
both oblate-like and prolate-like halos, Tp will be negative (positive) in the first and third (second 
and fourth) quadrants. This causes i to oscillate twice during one full precession of the LON (di/dt 
is positive [negative] in the fist and third [second and fourth] quadrants), i.e., the disk nutates as 
it precesses. 
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3.2. Torques by the Accreting Torus 

The accreting torus is inclined at a small angle (15°) relative to the halo equatorial plane, thus 
the torus has a negative J20, i.e., the most significant non-spherical component of the torus is the 
oblateness. Furthermore, since the mass distribution of the torus has a sinusoidal deviation from 
the equatorial plane with an azimuthal period of 27T, the torus has a non-zero J21 component (see 
Table [2] for these values at r = 2 and 4) , which is not present in triaxial halos and is much more 
important than the J22 component of the torus. 

Equations (IB8I) and (IB9I) show that di/dt and dX/dt are sinusoidal functions of A and imply 
that i and A will oscillate around certain equilibrium values. The circular velocity of the disk is 
nearly constant in the outer disk region, thus equations (IB8|h (IB9I) . and (IB2I) through (IB4I) imply 
that di/dt and dX/dt by the torus are approximately proportional to r, whereas di/dt and dX/dt 
by the halo decrease as r increases. 

3.3. Torques by the Inner Disk 

We adopt a simple two-disk (inner and outer disks) model to describe the effect of the self- 
gravity within the disk. Equations f)C3D and f)C2D imply that if A of the outer disk is ahead (in the 
sense of disk spin) of that of the inner disk (i.e., a leading spiral of the LON), di/dt is positive, and 
vice versa. Equations (IC4I) . (IC1I) and (IC2I) state that dX/dt by the inner disk is negative on average 
when the inclination of the inner disk, z^, is smaller than ~ 55°, and has a large A dependence 
when i of the outer disk is small. 



4. INITIALLY TILTED DISKS 

Here we present the results obtained from our simulations for the initially tilted disks and 
analyze them by comparing the axisymmetric and triaxial cases. 

4.1. Oblate vs. Oblate-like Halos 

Figure [2] shows the evolution of the initially tilted disk in an oblate halo (model TD-O). This 
simulation is set up in the same way as the oblate model of Ideta et al. (2000) in order to check the 
agreement between the two simulation results and to describe the results of our triaxial models in 
comparison to the axisymmetric ones. We find that our model TD-0 gives almost the same result 
as the oblate model of Ideta et al. First we briefly discuss the evolution of the disk in model TD-0 
in terms of torques following the analysis by Ideta et al. (2000). 

Oblate halos exert negative Tl on the disk, and the oblate halo adopted here has the absolute 
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precession rate that decreases as r increases. Thus the LON undergoes a retrograde precession 
with a shape of a leading spiral. Since the LON of the outer disk is ahead (in the sense of disk 
spin) of the LON of the inner disk, the inclination of the outer disk increases forming a type I 
warp (terminology used by Sparke & Casertano 1988), which bends upward away from the halo 
equatorial plane. This is because the inner disk behind the outer disk exerts a negative Tp on to 
the outer disk, which results in an increase of i (see §3.31) . However, a larger inclination causes a 
slower precession (\d\/dt\ oc cosz; see eq. |A17j ). thus the LON of the outer disk leads even further 
and does not have a chance to catch up the inner disk. As the LONs of the inner and outer disks 
become significantly misaligned, the disk is not observed as a warp anymore. 

This can be clearly seen in the i re \-\ diagrams (so-called "Briggs diagrams", Briggs 1990) 
in panel 6, which effectively shows the radial profiles of A and the inclination relative to that at 
the central region of the disk, i re \. Each circle represents an annulus of OARj between lRj and 
4.6^. The disk particles in some of the outermost annuli tend to spread out in the projected sky 
making them not recognized as a part of the disk in actual observations, so we use open circles to 
indicate the annuli that have a peak projected density less than 20 particles per 0.1 2 area. In such 
diagrams, an observable warp would be manifested as an extended alignment of filled circles. The 
filled circles are extended out to ~ 10° only until T ~ 200 and remain inside i re \ ~ 7° afterwards. 
Panel a indeed shows that a warp phenomenon is significant only until this epoch. 

Figure [3] plots the evolution of a disk in our moderately oblate-like halo, model TD-OL7 (note 
that all of our oblate and oblate-like halos have similar [a 2 + b 2 ]/c 2 ratios so that they have similar 
J20 values and thus similar precession rates). It shows that the triaxiality of the halo does not 
greatly change the global aspect of the evolution. The most significant difference between the 
oblate and oblate-like cases is that i and A oscillate periodically in the latter, which is the result 
of P- and L-component torques from the triaxial component of the halo, respectively — Tp and 
Tl are sinusoidal functions when J22 has a non-zero value (see eqs. |A14j and [A15j ). The phase 
of the oscillation coincides with the quadrant of the LON, which indicates that the oscillation is 
indeed caused by the triaxial component of the halo. As discussed in § [HI when A is in the first or 
third quadrant, the disk experiences a negative P-component torque and the inclination increases. 
Likewise, the inclination decreases when A is in the second or fourth quadrant. 

A similar evolutionary picture is observed even in the more extremely oblate-like case, model 
TD-OL13. Figure [4] shows that the global evolution of this model is very similar to that of model 
TD-OL7. 

In both models TD-OL7 and TD-OL13, the magnitude of the warp slightly oscillates, but the 
average magnitude of the warp is similar to that of model TD-O. To summarize, when J20 of the 
halo is negative, the warp in an initially tilted disk scenario is very weak except at the beginning, 
and even the triaxiality of the halo is not effective in increasing the magnitude of the warp. 
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4.2. Prolate vs. Prolate-like Halos 

The evolution of the disk in a prolate halo, model TD-P, is shown in Figure El Prolate halos 
exert positive Tl on the disk, and our prolate halos have an absolute precession rate that decreases 
as R increases. Thus the LON undergoes a prograde precession with a shape of a trailing spiral. 
Because the LON of the outer disk is behind (in the sense of disk spin) that of the inner disk, 
Tp exerted by the inner disk on to the outer disk is positive, and the inclination of the outer disk 
decreases (see §3.31) forming a type II warp, which bends down toward the halo equatorial plane. 
Since a smaller inclination causes a faster precession (\d\/dt\ oc cosz; see eq. |A17| ). the LON of 
the outer disk catches up that of the inner disk. If the LON of the outer disk outruns that of the 
inner disk, the outer disk will have a larger inclination and its precession rate will become smaller. 
Thus the self-gravity between the inner and outer disks plays a role of restoring force that leads 
to an alignment of the LONs between the inner and outer disks. The simulation shows that after 
about one revolution of the precession, the inner and outer disks settle into a stable equilibrium, 
and the LONs of the two disks remain close to each other afterwards, making the disk observed as 
a warp. 

This analysis can be generally applied to the evolution of the disk in a moderately prolate- 
like halo, model TD-PL7 (Fig. [6j). Note that all of our prolate and prolate-like halos have similar 
(a 2 + 6 2 )/c 2 ratios so that they have similar precession rates. Again, the most considerable difference 
from the prolate case is the oscillation in both i and A. Since the LON of the inner disk advances 
faster than that of the outer disk in the beginning of the simulation, the oscillation phase of the 
inner disk is ahead of that of the outer disk afterwards. Because of this phase lag, the inclination gap 
between the two disks increases and decreases periodically (without the phase lag, the inclinations 
of the two disks would oscillate with an almost constant gap). Therefore, as seen in panels a 
and 6, the magnitude of the warp oscillates with a period of precession as well. The warp nearly 
disappears at T ~ 250 and 400, and grows back to its peak at T ~ 350 and 500. Note that when 
the magnitude of the warp reaches its peak, it becomes larger than that of the prolate case with 
a similar prolateness (c 2 /[a 2 + b 2 ]). The i re \ value of the outermost annulus remains inside 15° in 
model TD-P, but it reaches near 20° at T = 350 and 500 in model TD-PL7. Thus the triaxiality 
of the prolate-like halo periodically suppresses and augments the magnitude of the warp. 

In a more extremely prolate-like halo (model TD-PL13), the disk evolves quite differently from 
the above two cases. As seen in panel c of Figure El the inclinations of both inner and outer disks 
increase fairly quickly while oscillating. This appears to be because of the sinz factor in di/dt term 
(eq. |A16| )-When the the magnitude of the triaxiality (J22) is large, the peak of the oscillating i 
can reach higher, and at a higher z, di/dt is larger. Thus when J22 is larger than a certain critical 
value, the peak i value in each oscillation gradually increases and eventually i becomes larger than 
90°. This happens at T ~ 500 for model TD-PL13, after which the LON evolves in the reverse 
direction because the cosi term in dX/dt (eq. |A17j ) becomes negative. 

The large J22 value and thus the large oscillation in i of model TD-PL13 causes the outer 
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disk to disperse so severely that it is not recognized as a part of the disk anymore after T ~ 450, 
making the size of the whole disk appear to be much smaller than in the beginning. The debris 
from the outer disk spreads out randomly around the inner disk after this epoch, and the warped 
configuration nearly disappears. 

Although panel c of Figure [7| shows the evolution up to T = 650, we have evolved model TD- 
PL13 up to T = 800 and find that after T ~ 400, i and A of the inner disk keep oscillating around 
90° and 360°, respectively, indicating that the disk plane tumbles about the minor axis of the halo. 
Such a phenomenon has been also observed in numerical and analytical studies on the dynamics 
of stars and gas in triaxial potentials (Aarseth & Binney 1978; Steiman-Cameron & Durisen 1984; 
Arnaboldi & Sparke 1994, among others): the particles in triaxial halos prefer the orbits whose 
principal plane tumbles about its minor axis. 

To summarize, when J20 of the halo is positive, an initially tilted disk can excite and sustain 
the warped configuration for a considerable time as far as the triaxial component of the halo is not 
too large, and a mild triaxiality can augment the peak magnitude of the warp. 

5. TILTED COSMIC INFALL 

Now we present and analyze our simulations for the cosmic infall scenario. In these simulations, 
the disk is initially in the equatorial plane of the halo and the dark matter torus, whose mass 
gradually increases over T = 100, is tilted at 15° relative to the halo equatorial plane such that its 
ascending LON coincides with the negative y-axis. 

5.1. Spherical Halo 

Figure [8] shows the evolution of the disk in a spherical halo (model CI-S). This model is to be 
compared with the axisymmetric and triaxial halo cases. As the halo in our model CI-S is perfectly 
spherical, the evolution of the outer disk is driven by the interaction with the torus and with the 
inner disk only. The tilted torus gives an effect similar to an oblate halo to the disk — it exerts a 
negative Tl on the disk and makes the disk precess in a retrograde sense about the axis normal to 
the torus plane. But unlike in model TD-O, the precession rate caused by the torus is an increasing 
function of R (see the Appendix A of Shen & Sellwood 2006 and our Appendix B, but note that 
the precession rate given in the former is along the torus plane while that in the latter is along the 
equatorial plane), resulting in a trailing spiral of LON. Thus the inner halo exerts a positive Tp on 
the outer disk, and the inclination of the outer disk decreases forming a Type II warp. 

The inclination relative to the equatorial plane evolves significantly with time (see panel c), 
but the inclination of the disk relative to the torus plane stays constantly (this phenomenon is to 
be compared with those in non-spherical halo models below). The increase and decrease of i are 
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well explained by equation (IB8I) . 

5.2. Oblate and Oblate-Like Halos 

Figure [9] shows the evolution of the disk in an oblate halo (model CI-O). Since the oblate halo 
does not have a J22 component, the change of i is caused only by the accreting torus. Once the i 
value departs from zero, A first moves toward —90° and then both i and A oscillate. The oscillatory 
behaviors are found to be drived by combined torques from the non-spherical halo and the tilted 
torus. Figure [10] shows four different sign combinations of di/dt and dX/dt from the combined 
torques, and the arrangement of the sign combinations is in such a way that the evolution of a 
point on the i-X plane follows a round, counterclockwise revolution about the equilibrium point 
(the point where both di/dt and dX/dt are equal to zero). Note that the center of A oscillation, 
—90°, coincides with the longitude of the ascending node of the torus, but the inclinations of both 
outer and inner disks are somewhat smaller than that of the torus. 

Inclinations of both inner and outer disks oscillate, but since the equilibrium point is located at 
a higher i for a larger r, the mean i value of the outer disk is higher than that of the inner disk, and 
thus the disk forms a type-I warp, the opposite of the spherical halo case. Thus in oblate halos, the 
magnitude of the warp or i re \ is determined by the radial gradient of the equilibrium inclination, 
and is found to be smaller than that in the spherical halo model. Unlike in model TD-O, the A 
values of the inner and outer disks in model CI-0 are not significantly separated. This is because 
in model CI-O, the inclinations of both inner and outer disks are so small that the cosz factors in 
equations (IA17I) and (IB9I) are not so effective in making the two A's separated as in model TD-0 
(i.e., the cosz values of the inner and outer disks are similar). 

Figure [TT] shows the evolution of the disk in one of our oblate-like halos (model CI-OL7). 
In this model, the shorter of the two halo axes in the equatorial plane aligns with the LON of 
the torus, thus the oblate-like halo induces a positive di/dt when A is in the fourth quadrant 
(270° < A < 360°) and a negative di/dt when A is in the third quadrant (180° < A < 270°). Such A 
dependence of di/dt is opposite to that by the torus, but the contribution on di/dt by the oblate-like 
halo is much smaller than that by the torus because of the small value of sinz in equation (IA16I) . 
The average inclinations of both inner and outer disks are smaller than those in the oblate halo 
model because the equilibrium inclination is smaller for the oblate-like halo model ( J22 cos 2A in eq. 
|A17j is negative when A is near the equilibrium A, 270°; see Fig. [12]). Other than these lowered 
equilibrium inclinations for both inner and outer disks, the overall evolution of the disk in model 
CI-OL7 is very similar to that in model CI-O. 

We find that the evolution of the disk in model CI-OL7b (figure not presented), where the 
longer, not the shorter, of the two halo axes in the equatorial plane aligns with the LON of the 
torus, is very similar to that of model CI-O as well, implying that the orientation of the tilt of the 
torus relative to the orientation of the oblate-like halo is not an important factor in determining 
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the configuration and the magnitude of the warp. The only noticeable difference between models 
CI-OL7b and CI-OL7 is the higher average i values for both inner and outer disks of the former 
model because of the positive J22 cos 2A value in equation fjA17D for A ~ 270°. Note that, however, 
all of our oblate and oblate-like halo models have nearly the same magnitudes of the warp (i.e., 
nearly the same i re i). 

5.3. Prolate and Prolate-Like Halos 

Figure fT3l shows the evolution of the disk in a prolate halo (model CI-P), which is rather 
different from those in spherical and oblate halos. In a prolate halo, the magnitude of the warp is 
much larger, and the outermost part of the disk is even detached from the rest of the disk when 
the LONs of the outermost part and the main disk are apart by ~ 180°. 

This rather large difference is thought to be because of the strong radial dependence of the 
combined dX/dt from the halo and the torus. Figure fT4l shows that the combined dX/dt for the 
prolate halo has a significantly stronger radial dependence than that for the oblate halo. Thus unlike 
in the oblate halo, the outer disk in the prolate halo can evolve separately from the inner disk from 
the beginning of the simulation. The LON of the inner disk oscillates around the equilibrium A, 
90° (see Fig. [T5]h whereas the LON of the outer disk, which is more affected by the self-gravity 
within the disk than the inner disk is, appears to experience more complicated evolution. In the 
beginning of the simulation (T < 150), the combined dX/dt on the outer disk from the halo, torus, 
and the inner disk is nearly zero, so A of the outer disk barely evolves. But once A's of the inner 
and outer disks are apart by more than 90°, the second term in equation (IC4I) becomes negative. 
This then breaks the subtle equilibirum and makes the outer disk precess in a retrograde fashion. 

Unlike in model TD-P, the A values of the inner and outer disks in model CI-P are well 
separated. In model TD-P, the cosi factor in dX/dt (eq. |A17j ) causes the outer disk, which has a 
smaller i than the inner disk, to catch up the A of the inner disk. But in model CI-P, the inclinations 
of the inner and outer disks are both so small in the beginning of the simulation that the cos i factor 
in the combined dX/dt from the halo and the torus is not so effective in making the two A's move 
together (|dcos i/di\ is smaller for small z's). 

Because A of the inner disk is initially ~ 90° ahead (in the sense of disk spin) of that of the 
outer disk, the inner disk experiences a positive di/dt from the outer disk (see eq. |C3| ). thus has 
an equilibrium i value that is higher than the one determined by considering the torques from the 
halo and torus only as in Figure [T5l The i value of the outer disk increases in the first half of the 
simulation and decreases in the second half of the simulation because the combined di/dt by the 
halo and the torus is positive when A is in the fourth quadrant and negative when A is in the third 
quadrant (see Fig. [T5l) . 

While the evolution of both inner and outer disks in spherical and oblate halos can be explained 
by the combined torques from the halo and the torus only, in a prolate halo, the self-gravity within 
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the disk plays a non-negligible role in the i and A evolutions of the disk. This is because in a prolate 
halo, 1) the di/dt term from the inner disk is not negligible as the LONs of the inner and outer 
disks are well separated from each other (see eq. |C3| ), and 2) the dX/dt term from the inner disk 
is relatively more important as the J20 values of the halo and the torus have opposite signs. 

Figures [16] shows the evolution of the disk in one of our prolate-like halo models (model CI- 
PL7), where the shorter of the two halo axes in the equatorial plane aligns with the ascending 
LON of the torus. The evolution of CI-PL7 is generally similar to that of CI-P except that the 
magnitudes of the warp and the detachment of the outer disk less significant than those in model 
CI-P. This is because the J22 component of this triaxial halo induces a negative di/dt on to the 
outer disk when A is in the fourth quadrant (i.e., during the first half of the simulation), causing 
the peak i that the outer disk acquires to be smaller than that in model CI-P. 

Figures [17] shows the evolution of the disk in another prolate-like halo (models CI-PL7b), 
where the longer of the two halo axes in the equatorial plane aligns with the ascending LON of the 
torus. The second term in equation (IA17I) now causes a positive (instead of negative as in model 
CI-PL7) dX/dt for the outer disk when A ~ 270°, and this apparently balances the negative dX/dt 
from the inner disk when the A's of the inner and outer disks are apart by more than 90° so that A 
of the outer disk stays near 290° after T ~ 200. At this equilibrium A, the combined di/dt on the 
outer disk from the halo and the torus is positive (see Fig. [T8|) . Consequently, i keeps increasing 
even after the outer disk is completely detached from the inner disk, and eventually the system 
resembles polar-ring galaxies. 

While the evolutions of the disks in oblate and oblate-like halos are alike and can be described 
rather simply and coherently without a consideration of the self-gravity within the disk, those in 
prolate and prolate-like halos sensitively depend on the relative importances between the torques 
from the halo, the torus, and the disk itself, as well as on the relative orientation between the halo 
and the torus. The two-disk analysis that we adopted to describe the effect of the self-gravity is 
found to give qualitatively correct descriptions, but does not give quantitatively accurate estimates 
for the torques. Thus it is rather difficult to generally describe the evolution of the disk when the 
self-gravity is not negligible, which is the case of the disks in prolate and prolate-like halos. 

6. SUMMARY AND DISCUSSION 

We have performed a series of simulations of galactic warps for two scenarios: 1) stellar disks 
that are initially tilted relative to the equatorial plane of triaxial halos, and 2) stellar disks in 
axisymmetric or triaxial halos with late infall of dark matter that is modeled as a torus tilted 
relative to the equatorial plane of the disk and the halo. 

In our initially tilted disk (TD) simulations, the triaxiality in the oblate halo does not signif- 
icantly alter the characteristics and magnitude of the warp, although the latter oscillates as the 
disk precesses and is smaller for a larger triaxiality. Moderate (7 %) triaxiality of the prolate halo 
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results in the warp similar to that in the axisymmetric prolate halo with an oscillatory behavior, 
while a larger (13 %) triaxiality makes the disk precess about the minor axis of the halo, instead 
of the major axis that the disk initially precessed about, and considerably reduces the magnitude 
of the warp. 

In our cosmic infall (CI) simulations, the magnitude of the warp in the non-spherical halo 
is determined by the difference in the "equilibrium z" between the inner and outer disks. While 
the magnitude of the warp in oblate and oblate-like halos is smaller than that in the spherical 
halo, in prolate and prolate-like halos the warp can grow larger and the outermost part of the 
disk can even detach from the main disk body. For a prolate-like halo system where the ascending 
LON of the dark matter torus aligns with the longer of the two halo axes in the equatorial plane, 
the detachment becomes more extreme and the system resembles polar-ring galaxies. Note that, 
however, cosmological simulations show that the angular momentum vector of the halo preferentially 
aligns with the minor axis of the halo (see the references below), thus the disk morphologies found 
in our prolate and prolate-like halo simulations (detached outer disk and polar-ring configuration) 
are not expected to take place often. 

Regardless of the warp scenario, the triaxiality of the halo causes i and A to oscillate (i.e., it 
causes the disk to nutate), the differential nutation between the inner and outer disks generally 
makes the magnitude of the warp slightly diminish and fluctuate, but this could also cause the 
outermost part of the disk to completely detach from the main body if this effect is boosted by the 
torque from the dark matter torus as in model CI-PL7b. 

In the Appendices we derived the approximate formulae for the torques exerted on the disk by 
the triaxial halo and the dark matter torus, and with these formulae, we have successfully described 
the i and A evolutions of the disks. The characteristics of the halo and the torus are represented 
by three parameters J20, J2I5 an d J22, which can be useful in predicting the evolution of a disk in 
various different halo and torus models, or in estimating the relative importance between the halo 
and the torus. 

The techniques used in deriving these formulae could be applied for halos with more com- 
plex structures. Bailin & Steinmetz (2005) found from their ACDM Af-body simulations that the 
principal axes of the galaxy dark matter halos are internally misaligned (i.e., the directions of the 
principal axes are functions of galactocentric radius). In such systems, although the stellar disk will 
quickly be in equilibrium with the inner halo, it may experience torques from the misaligned outer 
halo for an extended period, and the higher order terms than J21 and J22 might be important. 
Such configuration could be responsible for U- or L-shaped warps, and for these cases, each radial 
bin of the disk could be modeled with two broken half-rings, instead of a whole ring. 

The main goal of the present study is to understand the fundamental dynamics between the 
disk and the non-spherical halo (with and without the dark matter torus) , thus we have considered 
models that result in relatively prominent warps (a large initial tilt angle of the disk and a large 
non-sphericity of the halo). However, most of the observed edge-on disk galaxies have only a hint 
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of warped structure, if any, near — 4-5 in their stellar disks. For example, Ann Sz Park (2006) 
found that out of their 325 sample galaxies 165 (51 %) have an S-shaped warp, and 142 (86 %) of 
these S-shaped warps have warp angles a smaller than 5° (a is the inclination angle of the tip of 
the outermost isophotes from the mean disk major axis, and is very close to our i re i). 

Our simulations for the initially tilted disk scenario show that prolate and prolate-like halo 
models can excite more prominent warps than oblate and oblate-like halo models. However, this 
does not necessarily imply that, as Ideta et al. (2000) interpreted, oblate and oblate-like models 
are not appropriate to explain the observed warps-they may not be able to sustain a prominent 
warp, but surely can resemble the majority of the observed warps, i.e., those with small a values. 
On the other hand, prolate and prolate-like models can explain warps with both large and small 
a values by having different initial tilt angles or halo axis ratios. Figure [19] shows that i re \ can 
become smaller if prolate and prolate-like models have an initial tilt angle of 10°, instead of 30°, 
or a halo axis ratio of y/ab/c = 0.85, instead of 0.75. 

Our simulations for the tilted cosmic infall scenario also show that prolate and prolate-like halo 
models can excite more prominent warps than oblate and oblate-like halo models. The outermost 
disk in those models can detach from the main disk body or even form a polar ring at least for some 
time, but such systems are barely observed. The magnitude of the warp is suppressed in oblate 
and oblate-like halos, and is even smaller than in the spherical halo. Thus for the cosmic infall 
scenario, the way to obtain a nice and clean, prominent warp would be to have a (nearly) spherical 
halo and/or dark matter infall that gives a large J21 value. 

These results appear to be consistent with the recent A cold dark matter cosmological simu- 
lations. Bailin et al. (2005) found in their hydrodynamic simulations that the (inner) halo minor 
axes are nearly aligned with the disk axis, and Bett et al. (2007) analyzed the Millennium sim- 
ulation (Springel et al. 2005) to find that the majority of halos have their angular momentum 
vector aligned with their minor axis. Both of these, which confirm the findings of earlier studies by 
Dubinski (1992) and Warren et al. (1992), imply that most of the disks are situated in an oblate 
or oblate-like halo, and can explain why 1) most of the observed warps are quite weak (warps are 
suppressed in oblate and oblate-like halos) and 2) disks with detached outer part or polar-ring 
systems are not observed or rare. 

We are grateful to Juhan Kim, Woong-Tae Kim, Jounghun Lee, Juntai Shen, and Linda 
Sparke for helpful discussion. We thank the anonymous referee for valuable comments. M. J. 
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for the Structure and Evolution of the Cosmos (ARCSEC) of the Korea Science and Engineering 
Foundation through the Science Research Center (SRC) program. Simulations presented in this 
paper were performed on the GRAPE cluster at Kyung Hee University and on the Linux cluster at 
Korea Astronomy & Space Science Institute (KASI), which was built with the fund from ARCSEC 
and KASI. 
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A. Approximate Torques by the Triaxial Halo 

Here we derive approximate torques exerted on a ring of disk material by triaxial halos following 
the procedure in Steiman-Cameron & Durisen (1984), who considered the halo mass distribution 
inside the ring only We extend their work and consider the halo mass distribution both inside and 
outside the radius of the ring. We then derive the rates of changes in inclination and LON of the 
ring due to the non-spherical components of the halo. 

Gravitational potential from a mass distribution p can be written as an expansion in spherical 
harmonics Yj m , 



(Al) 



i= 0m =-i M + l Lr Jo Jr r 

Here pi m are the ^-coefficients of the shell lying between r and r + Sr, 

Plmir') = I Y^O, <f>)p(r>, 0, <f>)dto, (A2) 

where Yf m are the spherical Legendre conjugate functions, and d£l is the solid angle element. Then 
the first-order corrections to the potential can be written as 

G 

$i = -- [J 20 (3cos 2 e - 1) + 3J 22 sin 2 9 cos 2(f\ , (A3) 



where J20 and J22 are defined by 



J2 (r) = J$(r) + r 2 t(r) 

Mr) = J%\r) + r 2 ?{r) (A4) 



and 



1 



r =r 



4o(r) = -W (2z' 2 -x' 2 -y' 2 )pdV (A5a) 

T Jrr'—Pi 



rr =00 

JZo\r) = r 2 {2 z' 2 -x' 2 -y' 2 )-^dV (A5b) 

J r'=r f 



1 



r —r 



4f(r) = -3 / (x' 2 -y> 2 )pdV (A6a) 

rr' =00 

J2?(r) = r 2 ^-y'^P- d v'. (A6b) 

J r' =r f 

Oblate halos have negative J20 values while prolate halos have positive J2o's. Also note that oblate 
and prolate halos have J22 = 0, and triaxial halos with a > b (b > a) have positive (negative) J22's. 
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Now with equations (A16) to (A18) of Steiman-Cameron & Durisen (1984), the average first- 
order potential corrections for a ring of disk material that has a radius of r, inclination z, and the 
longitude of ascending LON A become 

G 

($i> = -- [J 2 o(3sin 2 i - 2) + 3J 22 cos 2A sin i] . (A7) 
8 

If the spin of the ring is much faster than the precession or the inclination change of the ring, the 
Lagrange equations of motion for the ring become 

n _ i a<*i> 



dt rV c sin i dX 

dX _ 1 <9($i) 
dt rV r sin i di 



(A9) 



(Habe & Ikeuchi 1985), where V c , the circular (spin) velocity of the ring, is defined to be positive 
when the ring spins counterclockwise when viewed from the positive z-axis, and the Hamiltonian 
equations of motion for the ring yield 

= ^ , 

oi 

h - -?^y> (ah) 

(Arnaboldi & Sparke 1994). Since the L-axis coincides with the LON, we have 

T L =Pi (A12) 

and since the P-component angular momentum of the ring is — mrV c smi, we have 

di 

T P -mrV c cos i— (A13) 

(see $3] for the definitions of L- and P- axes). Now the L and P components of the torque become 

0(<fri(M,A)) 

TL = — di — 

3G 

= ( J20 + J22 cos 2A) sin % cos i ( A14) 

. a(^(r,z,A)) 

Tp -cot 2 — 

aA 

= J 2 2sin2Asin2z. (A15) 

8 

Finally, from equations (|A8I) , (|A9I) , (|A14I) , and ()A15I) , the rate of inclination change and the 
precession rate can be given by 

di Tp 



dt rV c cos i 
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dX 
dt 



3G_ 

4rV c 



J22 sin 2A sin i 



rV c sin i 
3G 



4rV r 



( J20 + J22 cos 2 A) cos i. 



(Me) 



(A17) 



B. Approximate Torques by the Dark Matter Torus 

Here we derive approximate torques exerted on a ring of disk material by the accreting torus 
of a mass M t that is inclined at i t from the equatorial (x-y) plane of the halo with the ascending 
LON aligned with the negative y-axis. We assume that the torus has an inflnitesimally small cross 
section and a radius of Rt that is larger than that of the ring. 

The potential by the torus does not have / = 1 terms as in the case of triaxial halos, but it 
does have I = 2, m = ±1 terms besides I = 2, rn = 0, ±2 terms. Then the first-order corrections to 
the potential can be written as 

G 

$1 = -- [J 2 o(3cos 2 6 - 1) + 6J21 sin 29 cos (j) + 3 J22 sin 2 9 cos 20] . (Bl) 

We find that the coefficients J20, J2I5 an d J22 have the following forms for the torus: 

M-l r 2 

J2 ° (r) = ^^3( 3sin2 ^- 2 ) ( B2 ) 

Mt r 2 

Jn(r) = Xif sin2 ^ (B3) 

M, r 2 

= -^-^sin 2 *,. (B4) 

Now with equations (A16) to (A18) of Steiman-Cameron & Durisen (1984), the average first- 
order potential corrections for a ring of disk material that has a radius of r, inclination z, and the 
longitude of ascending LON A become 

G 

($1) = [j 2 o(3sin 2 i - 2) - 6 J 2 i sin A sin 2i + 3 J22 cos 2 A sin 2 i] . (B5) 

8 

Then from equations (IA14I) and (IA15I) . the L and P component torques are given by 

3G 

Tl = — ( J20 sin 2i — 4 J21 sin A cos 2i + J22 cos 2A sin 2i) (B6) 
8 

3G 

T P = — (4J 2 i cos Acos 2 z + J 2 2sin2Asin2z), (B7) 
8 

and as in Appendix [Aj the rate of inclination change and the precession rate are given by 
di 3G 



dt 4rV c 



(2 J21 cos A cos i + J22 sin 2A sin i) 
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3G 

2 J21 cos A cos i (B8) 



4rV, 



C 



d\ 3G ( . A cos 2i 

J20 cos i — 2 J 21 sin A— — - + J 22 cos 2 A cos 1 



dt ArV c \ sinz 

3G ( e cos2i\ 

J20 cos z — 2 J 21 sin A— — r • (B9) 



4rl/ c \ sinz 

Note that for our torus, J22 terms are negligible compared to the other terms (see Table [2]). 

Shen & Sellwood (2006) derives the precession rate of a disk due to an accreting torus by 
assuming that the torus lies in the equatorial plane and the disk is inclined at a certain angle 
relative to the equatorial plane (note that our i and A are defined to be the inclination and the 
LON respectively in the coordinates where the torus is inclined relative to the equatorial plane; 
when the torus lies in the equatorial plane, the inclination of the disk does not change and only A 
of the disk evolves). The rates of changes in i and A can be also derived from the precession rate 
by Shen & Sellwood, 

d\ f 3GM t r 



dt 4R*V C 



cosz', (BIO) 



where X is the LON of the disk along the torus plane and i' is the inclination of the disk relative 
to the torus plane. By transforming the coordinates from the torus plane to the equatorial plane, 
one obtains the following rates of changes: 

di . d\ f 

— — cos Asm it—— (Bll) 
dt dt y J 

d\ , v dW , v 

— = (cos it + cot i sin A sin i t ) — — . (B12) 

CLL CLL 

The cosi f term in d\ f /dt can be transformed to the equatorial coordinates by 

cos i f — — sin i sin A sin i t + cos i cos i t . (B13) 



C. Approximate Torques by the Inner Disk 

Here we estimate approximate torques on the outer disk by the inner disk using the results 
obtained in Appendix [Bl If the inner disk is considered in place of the torus, equations (jB2l) and 
(1B3P can be written as 

J20W = ^%W^- 2 ) (CI) 

Zi IT 

J 2l(r ) = M^tik sln2l . di (C2 ) 
4 r a 



and the rates of changes become 

di 3G 



dt 4rV c 



2 J21 sin(A — A^) cos i (C3) 
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dX 3G ( _ x cos2i\ 

J 2 o cos z + 2 J 2 i cos(A - X id ) — — r , (C4) 



4r y c \ sin i 

where the variables with a subscript id are for the inner disk and those with no subscripts are for 
the outer disk. 
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Table 1. Simulation Parameters 
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7.5 
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Ma 

IN (J 
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Oblate-like, 7 % 
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TD-OL13 


11.3 


8 7 
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Oblate-like, 13 % 


No 

1NU 


TD-P 


7.5 




10 


Prolate 


Ma 

IN (J 


TD-PL7 


8.03 


u.yo 


10 


Prolate-like, 7 % 


No 

1NU 


TD-PL13 


8.48 


6.53 


10 


Prolate-like, 13 % 


No 


Tilted Cosmic 


Infall 










CI-S 


10 


10 


10 
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Table 2. J Values of Our Models 



Model 


r = 2 




r = 4 




^20 ^21 




^20 ^21 


^22 


Halo 










TD/CI-0 


-5.05E-2 




-6.67E-2 




TD/CI-OL7 


-5.00E-2 


1.22E-2 


-6.62E-2 


1.62E-2 


TD/CI-OL13 


-4.90E-2 


2.28E-2 


-6.49E-2 


3.09E-2 


TD/CI-P 


5.84E-2 




7.61E-2 




TD/CI-PL7 


5.89E-2 


1.43E-2 


7.66E-2 


1.81E-2 


TD/CI-PL13 


6.01E-2 


2.66E-2 


7.78E-2 


3.36E-2 


Torus 










All 


-7.31E-3 1.02E-3 


-2.72E-4 


-2.93E-2 4.07E-3 


-1.09E-3 
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Fig. 1. — Definitions of L- and P- axes. The positive L-axis coincides with the line of ascending 
node (LON) of the ring, and the P-axis is normal to the L-axis in the halo equatorial plane such 
that L x P aligns with the positive z-axis. Note that both L- and P- axes always lie in the x-y 
plane. 
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Fig. 2. — Disk Evolution of model TD-O. a) Projected density distribution of the simulation 
particles as seen from the ascending line of node (LON) of the disk. The abscissa aligns with 
the x-y plane. The color bar shows the level of the surface density in units of # particles per 
projected area of 0.1 2 . b) i re \-\ plot (the 'Briggs diagram') for 9 equal-size radial bins between 
r = 1.0 and 4.6 (thus each bin has a width of 0.4). The distance from the plot center corresponds 
to i re l, the inclination relative to the central region of the disk, and the azimuthal angle from the 
positive x-axis corresponds to A, the longitude of the ascending LON. Each concentric circle marks 
10°. The inclination of the disk center is defined to be the average i of the two innermost radial 
bins(r = 1.0-1.8). Open symbols are for the radial bins whose peak projected density is smaller 
than 20 particles per 0.1 2 area, c) i evolution of two radial bins that represent the inner and outer 
disks (the radial ranges of these bins are indicated in the plot), d) A evolution of the same two 
radial bins as in panel c). 
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Fig. 3. — Same as Figure EJ but for model TD-OL7. 
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Fig. 4. — Same as Figure EJ but for model TD-OL13. 
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Fig. 5. — Same as Figure El but for model TD-P. 
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Fig. 6. — Same as Figure EJ but for model TD-PL7. 
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Fig. 7. — Same as Figure EJ but for model TD-PL13. 
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Fig. 8. — Same as Figure El but for model CI-S. The 10 equal-size radial bins for panel b) range 
from r = 1.0 to 6.0, each bin with a width of 0.5, instead of 0.4. 
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Fig. 9.— 



Same as Figure EJ but for model CI-O. 
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Fig. 10. — di/dt = (thick lines) and dX/dt = (thin lines) curves for the inner disk (r = 1.5; upper 
panel) and the outer disk (r = 5.0; lower panel) of model CI-O. For di/dt and d\/dt, contributions 
from both halo (eqs. [ATH] & [ATT] ) and torus (eqs. [B8j fc |B9j) are considered. Tickmarks 
are drawn on the side of negative values (negative time derivatives). The upward arrow indicates 
that the variable increases in that region, while the downward arrow indicates that the variable 
decreases. The dashed circle with an arrow is centered at the "equilibrium point" (see the text for 
its definition) and illustrates how (i, A) of the disk would circulate around the equilibrium point. 
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Fig. 11. — Same as Figure El but for model CI-OL7. 




Fig. 12.— 



Same as Figure [THl but for model CI-OL7. 
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Fig. 13. — Same as Figure El but for model CI-P. 
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Fig. 14. — Radial profiles of dX/dt for models CI-0 (solid line) and CI-P (dashed line) when i = 10°. 
Contributions from both halo (eq. [A17] ) and torus (eq. |B9j ) are considered. dX/dt of model CI-P 
varies more significantly than that of model CI-0 between r = 1 and 5, which is thought to be the 
main cause of the independent evolution of the outer disk from the inner disk. 
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Fig. 15. — Same as Figure [TU1 but for model CI-P. 
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Fig. 16. — Same as Figure El but for model CI-PL7. 
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Fig. 17. — Same as Figure El but for model CI-PL7b. 
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Fig. 18. — Same as Figure fTUl but for model CI-PL7b. 
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Fig. 19. — i evolution of the two radial bins for models that are variants from models TD-P7 (left 
panels) and TD-PL7 (right panels). The upper panels are the models with an initial disk inclination 
of 10° instead of 30°, and the lower panels are the models with an bulge axis ratio y/ab/c of 8.5 
instead of 7.5. These variant models result in smaller magnitudes of warping than our canonical 
models, TD-P7 and TD-PL7. 



